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I. 



INTRODUCTION 



Time-non-local (TNL) Lagrangians are generalized Lagrangians which contain integral terms 
involving the whole history of the system. These Lagrangians are of interest in different fields of 
theoretical physics, from string theory [lH3], to regularized QFT [11[5], semiclassical gravity [6], 
non-Markovian open quantum systems [7H9] , collapse models |10H13j , and recently light harvesting 
in photosynthesis process [TH [TB] . 

A major and typical difficulty one encounters when working with TNL Lagrangians is the 
following. Consider for example: 



where q and q are respectively the position and velocity of a particle, and a{t, s) a memory kernel [71- 
I12t [T^. Suppose one wants to compute the associated Hamiltonian, and needs to transform the 
dependence on (g, in a dependence on (or viceversa). The first guess would be to use the 

standard definition given by the Lagrange formalism: P = However, one can immediately see 
that, since in Eq. ([T]) the velocity q enters the integral term, the differentiation of L with respect 
to q does not make sense, and hence one does not know how to compute p. One therefore needs to 
develop a new formalism which allows for the analysis of TNL Lagrangians. In doing so, one has 
always to make sure that, in the limit a(t, s) = 5{t — s), i.e. when the Lagrangian becomes local, 
the standard results are recovered. 

Two are the main mathematical frameworks which have been developed, to deal with TNL 
Lagrangians: the Ostrogradsky formalism [16^ [T7j and the one proposed by Llosa and collabora- 
tors |18| [19] . The Ostrogradsky method is based on the assumption that some TNL Lagrangians 
can be represented as the limit of a sequence of finite higher derivative local Lagrangians. In this 
way one can apply an extension of the standard Euler-Lagrange formalism and write the equations 
of motion. This is reminiscent of the well known property in probability theory, according to which 
a non-Markovian process can be seen as the projection of a Markovian process defined on a higher 
dimensional space. It is worthwhile mentioning that the Ostrogradsky formalism cannot be applied 
to any TNL Action p^fTS]. 

Another interesting approach is the one proposed by Llosa and collaborators \18\ [T9] . They 
introduce in the Lagrangian a new parameter, in such a way that the evolution derived by this 
modified Lagrangian becomes local in the time variable, but non-local in the new parameter. In 
this way one can apply the standard Lagrange formalism with respect to the time variable, and 




(1) 
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then translate the results back to the non-local case. However, unlike the standard time-local case, 
the authors assume the Action to be a functional of the position only. Moreover, the formalism can 
be applied only to those Actions for which the functional and time derivatives commute [TU] . 
This approach allows for some reasonable mathematical treatment, but takes away from physical 
intuition. Other formalisms which have been proposed are less convincing or still debated |20| 121]. 

Aim of this Letter is to introduce a novel TNL Euler-Lagrange formalism which generalizes, 
in a somehow natural way, the conventional formalism, which is mathematically straightforward 
to handle. Moreover, such a formalism has the advantage of naturally emerging from the Hamil- 
ton's principle, and in the time-local limit automatically reproduces the standard Euler-Lagrange 
formalism. 



II. THE FORMALISM 



We consider a classical system described by the following TNL Lagrangian: 

L[q, q, t] = L^i^iq, q, t) + LtnlI^, q, t] , (2) 

where L^^i^iq, q, t) is the standard time local (TL) Lagrangian and LtnlIq, t] is the time non-local 
part, i.e. a functional of q and q which accounts for the whole past history of the system: 

^TNL [q,q,t]= f{q,q, s) ds , (3) 
Jo 

with f{q,q,s) a generic function. Eqs. ([2])-([3]) refer to a one-particle system; the generalization 
to many-particle systems is straightforward. Here and in the following, we denote functions with 
parenthesis and functionals with squared brackets. We require the Hamilton's principle to hold 
true, which means that the first variation 6S of the Action functional 

S[q,q] = / L[q,q,s]ds, (4) 
Jo 

vanishes. Note that in the standard case, that is when Ltnl[^, <?, = so that the Action ^ 
becomes local, one has 

SS[q, q]= n ^J^l4^6q + ^^4^13^] ds . (5) 
Jo \ dq dq J 

Integrating this expression by parts and requiring the variation to be zero, one obtains the well 
known Euler-Lagrange equations 1221 EH]. However, if LtnlIq, 7^ the standard Lagrange 
formalism cannot be directly applied, because the Action q] contains a double functional. 
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Therefore, differentiation with respect to q (or q) at a fixed time does not make sense anymore. 
This is the reason why one needs to go beyond the standard formalism. 
A way to circumvent this difficulty is to use the identity 

ft 



(6) 



where S/6q{s) denotes the functional derivative with respect to q{s). Integrating by parts Eq. ^ 
one obtains: 

Comparing this result with the local one, one then naturally defines the TNL Euler-Lagrange 
equations as follows: 

Some comments are at order. This equation of motion is an integro-differential equation which 
depends on the whole time interval [0,t]. As we will see in the following example, the correct 
interpretation of the Eq. ([s]) will require special care. Secondly, one can easily check that when the 
Lagrangian is local, i.e. LtnlI?, q, t] = 0, Eq. ^ reduces to the standard Euler-Lagrange equations. 
As a third comment, our formalism looks similar to the one used for (space) non-local fields in 
Quantum Field Theory. However, in that case fields are non-local only in the space variables and 
not in the time one j25l 126] . The fact that in our case non-locality affects the evolution variable, 
drastically changes the physical situation. We also stress that the presence of the Action S[q, q] 
instead of the Lagrangian in Eq. ([s]) is crucial, otherwise the formalism would be inconsistent (the 
functional derivative would generate Dirac deltas). 

We carry on the analogy with the standard formalism, introducing a generalized momentum: 

P(s) = ^|||i, V.G[0,t], (9) 

and defining the following generalized Hamiltonian functional: 

'H[q,p] = -S[q,q]+ f p{s)q{s) ds . (10) 
Jo 

The reason for introducing is that, as in the TNL Lagrangian formalism we need to express 

the functional equations of motion in terms of the Action instead of the Lagrangian, in the TNL 
Hamiltonian formalism we need to use the generalized Hamiltonian Hlqjp] instead of the standard 
Hamiltonian, otherwise, infinities would appear from the functional derivatives. We now evaluate 
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the first variation oiT-L[q,p], from which we will derive the equations of motion in the TNL Hamilton 



formalism. Using Eq. ( 10 ) one can write 



6n[q,p] = 1^ [-^-^j^^Qi^) - ^-^^Sq{s)+5pis)q{s)+p{s)6qis)^ ds 

= 1^ [-^-^^^l(') + ^Pi^M^)) ds , (11) 

where last line comes from Eq. ^ . Comparing this expression with the general form of the variation 
of n[q,p] [25j: 



one finds that 

dn[q,p] _ dS[q,q] Sn[q,p] _ 

6q{s) ~ 5qis) ' 6p{s) ^ 

This completes the TNL Hamiltonian formulation of the equations of motion. So far, we have 
developed the formalism under a rather formal point of view. We now show how this formalism 
applies to a specific, and very common in physics, example of TNL Lagrangian. 

III. APPLICATION OF THE FORMALISM 

We consider a classical particle of mass m evolving according to a general TNL second order 
Lagrangian, given by Eq. ^ with: 

LtUq, q, s) = "^qHs) + Aq{s)qis) + Bq\s) + Cg(s) + Dq{s) , 

LTNi.[q,q,s] = / a{s,r){Eq{s)q{r) + Fq{s)q{r) + Gq{s)q{r) + Hq{s)q{r))dr , (14) 
Jo 

where the coefficients A, . . . ,H are constants, and a(t, s) is a two-variable function. More general 
expressions for the Lagrangian can be considered, and the formalism still applies. First of all, we 
show that such a Lagrangian, for the purposes of our formalism, can be rewritten in a simpler but 
equivalent way. The reason is that, since the crucial quantity is the Action S, which is the integral 
of the Lagrangian, then Ltnl can be rewritten in different ways, without changing S. In particular, 
it is convenient to integrate by parts some of the terms containing double integrals. For example, 
the first term of Jq Ltnl^, q, s] can be rewritten as follows: 

/■* /■* / da(s,r) \ 

^ / qi^) / a(s,r)q(r) dr ds = E / q(s){a{s,s)q(s)—a(s,0)q(0)— / — ' — q(r) dr ] ds , 

Jo Jo Jo V Jo dr J 

(15) 
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and the first two terms on the right-hand-side can be absorbed in the term of L^l proportional to 
q{s). The same happens for the third and fourth terms of L^nl s]. Accordingly, without loss 
of generality, L[q, q, s] can be rewritten as follows: 

L[q, q, s] = "^iis) + A{s)q{s)q{s) + B{s)q\s) + C{s)q{s) + b{s)q{s) + F{s)q{s) a{s, r)q{r) dr , 

(16) 

where now the coefficients A{s) . . . F{s) are functions of time, which collect all the terms coming 
from the integrations by parts described here above. We can now use the formalism previously 
introduced, in particular Eq. ([s]) (or, equivalently, the generalized Hamiltonian formalism) to find 
the second order integro-differential equation of motion, which turns out to be: 

mq{s) = -A'{s)q{s) - C'{s) + 2B{s)q{s) + D{s) + F{s) [ a{s,r)q{r) dr + [ F{r)a{s,r)q{r) dr . 



Jo Js 

(17) 

This equation represents the generalization of the equation of motion of the standard local theory. 
The first thing one notices is that it displays a dependence on the future of s, not on its past, 
as one would naively expect by looking at the Lagrangian. The reason is that the problem was 
originally formulated with boundary conditions at the initial time and the final time t: finding 
the path that minimizes the Action, among all those which take the value go ^-t time and q at time 
t. While, in the standard time-local case, these boundary conditions can always be replaced with 
initial conditions on the position qQ and velocity go at time 0, in the TNL case this is not possible 
anymore, in general. This does not imply that the dynamics of the system necessarily depends on 
the future. It could, like in the Abraham-Lorentz equation [27] for a charged particle, which takes 
back-reaction into account, and displays what is called pre-acceleration; or like in the Weehler- 
Feynman formulation of classical electrodynamics [28]. But it could also not depend on the future, 
as it happens for non-Markovian stochastic processes, even if the dynamical equation formally 
does depend on the future. In this case, the equation does not correspond to the intuitive way of 
describing the dynamical evolution of the system. In general, there is no rule for discriminating 
which TNL equations can be re-written in such a way to show only a past-dependence, and which 
really contain a dependence on the future. From a mathematical (but also physical) point of view, 
this is the major difference between time- local and TNL dynamics, and why TNL dynamics are 
by far much more complicated. We will come back on this issue. As a final note, it is easy to 



check that in the limit F{s) = or a{t, s) — t- 6{t — s) Eq. (17) reduces to the standard equation of 
motion. 

In order to complete our general example, we compute the generalized Hamiltonian associated 
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to the TNL Action of Eq. (16). In particular, computing p{s) with Eq. ([9|, and substituting the 



corresponding expression, as weU as Eq. (16), in Eq. (10), one finds that: 
'H[q,p] 







^p\s) - ^qis)p{s) + - B{s)) q\s) - ^pis) 

2m m \ 2m J m 



+ 



A{s)C{s) - D{s] 



F{s)q{s) / a{s,r)q{r) dr 
Jo 



ds (18) 



m ) 2m 

Also in this case, one can easily check that in the time-local limit (a(s,r) — t- b{r — s)) this TNL 



Hamiltonian reduces to the standard one. Moreover, the transformation in Eq. (10) can be reversed 



and, substituting in it Eq. (13), one recovers S'[(7,(7] from Eq. (18) 



IV. EXAMPLE OF A TNL SYSTEM: THE TNL HARMONIC OSCILLATOR 



We now apply the TNL formalism to a specific physical situation in order to better understand 
the effects of the TNL terms on the dynamics of the system. In particular we consider a TNL 
harmonic oscillator, described by the following Lagrangian: 



-'^Qit) -kqit) ait,s)qis)ds., 



(19) 



where an integral term, which accounts for the memory effects, is added to the local Lagrangian. 
The equation of motion becomes: 

mq{s) = —kq{s) — k / a{s,r)q{r) dr . (20) 
Jo 

Again, we can see that the equation of motion formally depends on the future up to time t, other 
than depending on the past. To show that this is not necessarily a real physical dependence on 
the future, let us consider an exponential memory kernel 



a{t, s) = —e 



2 



(21) 



and show how the integro-differential equation can be transformed into an ordinary differential 
equation. (Note that by taking the limit 7—7-00 one recovers the time local case.) With this 



choice for the memory kernel, Eq. (20) becomes: 

mq{s) = —kq{s) 



kj 



--yli — s\ 



q{r) dr . 



(22) 



Note this is the same equation one obtains by applying the formalism of Llosa et al. [18] to the 



Lagrangian (19). Differentiating Eq. (22) with respect to s one finds 
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and differentiating it once more and using Eq. (22), one can show that q{s) satisfies the following 
fourth order differential equation: 



m 



+ 7^ 1 q{s) + —{k + k)q{s)=0. 
m 



(24) 



The first important message is that this is an ordinary differential equation, with no explicit 
dependence on the future (or past). This makes the point that, quite generally, one needs special 
care when interpreting the meaning of the solutions of integro-differential equations. 



The space of solutions of Eq. (24) is larger than that of Eq. (22). Therefore, one has to 



force suitable constraints on the solutions of Eq. (24) in order for them to solve Eq. (22) as well. 



In particular, since Eq. (24) is a fourth order differential equation, four boundary conditions are 



needed, to univocally determine the solution. Two of them are given by the boundary conditions of 
the variational problem, q{0) = q^ and q{t) = q. The dependence on the endpoint can be removed, 
like in the time local formalism, by replacing it with the initial condition on velocity g(0) = vq. The 



other two boundary conditions are given by the "consistence conditions" of Eq. (24) with Eq. (22), 



i.e. by asking the solutions of Eq. (24) to be solutions of Eq. (22) as well. In order to find these 



two constraints, one evaluates Eqs. (22)-(23) at time s = and then at time s = t, and compares 



all terms of these equations which are equal. The result is: 

g (0) - -m = 7 {m - -q{0) 
m \ m 

q{t) - -q{t) = -7 (q{t) - -q{t) 
m \ m 



(25) 
(26) 



The first consistency condition can be considered as an extra initial condition. The second con- 
sistency condition, instead, refers to the terminal time t. We are not aware of any way to replace 
it with an additional initial condition at time 0; accordingly, our original problem with boundary 
conditions at s = and s = t can only partly be reformulated as a problem with initial conditions 
at s = 0. But once again, we stress that this does not imply that the dynamics depends on the 
future. It is a consequence of the fact that the original problem was formulated with boundary 
conditions both at time and at time t. 



The solution of Eq. (24) reads: 



q{s) = ai^t sinhxis + a2^t coshxis + a^^^t sinha;2S + 04^4 cosha;2S , 



(27) 



where the subscript t denotes the dependence on the endpoint, due to the condition (26). The 



roots Xi and the coefficients ai^t are uniquely determined, and their explicit expression is given in 
the Appendix. We now discuss the physical behavior of the solutions. 



9 




Figure [T] shows the behavior of a TNL harmonic oscillator of mass m = 1, with initial position 
qo = and initial speed vq = 1. The strength of the coupling constants k, k, and the correlation 
function cutoff 7, have been chosen in such a way to highlight memory effects: k = 1, k = 10^, and 
7 = 1. We see that, as expected from a dynamics with memory effects, after an initial transient of 
few oscillations, the trajectory stabilizes and becomes harmonic. We notice that, if also the initial 
velocity is vq = 0, the oscillator stays at rest like in the time-local case. 

Figure [2] shows a comparison between the behavior of a standard harmonic oscillator (m = 1 , 

= 1) 7 = 00) and those of TNL harmonic oscillators (same m and vq) for different values of 
cutoff: 7 = 0.1, 0.3, 0.7. As we can see, for very shorts times, the behavior of the non-local 
oscillator follows that of a local one. Moreover, as expected, the bigger 7, the closer the non-local 
trajectories are to the local one. We further notice that the wavelength and the amplitude of the 
oscillations to which the TNL oscillator stabilizes are bigger than the corresponding local oscillator. 

We can then identify three main stages in the evolution of the TNL oscillator: 



1. Very-short time behavior: The TNL harmonic oscillator behaves as the local one. The reason 
is that the elapsed time is so short that there cannot be any memory of the past dynamics 



yet. In fact, in this stage the contribution coming from the integral term of Eq. (22) is 
negligible. 



2. 



Transient phase: In this stage, the oscillator behaves in an highly non harmonic way; this is 
when the TNL terms start affecting the dynamics. The length of this stage clearly depends 
on the non-local coupling constant k and memory kernel cutoff 7. 
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FIG. 2. Time evolution of q{s) of an oscillator with m = 1, for different values of cutoff 7 (dashed lines). A 
comparison with the standard harmonic oscillator (solid line) is given. The input data are: Xq = 0, vq — 1, 



3. Long-time behavior. The trajectory stabilizes to that of a local harmonic oscillator. The 
reason is that when the elapsed time is larger than the characteristic time of the memory 



kernel, the integral term of Eq. (22) reaches, for all practical purposes, an asymptotic stable 
value, and the subsequent dynamics looks "local-like". The long-time wavelength of the TNL 
oscillator is larger than that of a local one (7 = 00), because the effective force to which 
the oscillator is subject is weaker. The reason is that when q{s) is, like in the present case. 



an oscillating function, the value of the integral of Eq. ( 22 ) for any finite cutoff 7 is always 
smaller than the local case (7 — )• 00). 

These features of the TNL harmonic oscillator correspond to the physical intuition one would have 
about a standard harmonic oscillator interacting with a non-Markovian bath. 



V. CONCLUSIONS 



We have proposed a TNL formalism which allows for a general analysis of TNL Lagrangians, 
and we have derived the corresponding Hamilton formalism. As we have discussed, to study TNL 
dynamics starting from their Action, one has to solve a problem with boundary conditions at time 
and time t, which in general cannot be reduced to an initial-value problem (contrary to the standard 
time-local case, where this can always be done). This can be embarrassing in some cases, since 
the solution of the problem turns out to formally depend on the future. However, in important 
situations — like the path integral formalism — one needs only to solve a problem with boundary 
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conditions at the initial and final time, without any need to transform in into an initial-value 
problem. Therefore, the TNL formalism can be applied to such a formalism without difficulties, 
other than technical ones related to the evaluation of the integrals. 
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VI. APPENDIX 



We provide the explicit expressions for the roots Xi and the coefficients of Eq. (27), which 



determine the solution of Eq. ( 24 ) for the given initial conditions. The roots Xj are the solutions 



of the characteristic polynomial associated to Eq. (24) and they read {i = 1,2): 



2 



m 



k 
m 



+ T 



m 



(28) 



Since the coefficients have the same denominator d, is it useful to define four new coefficients 
bi as follows: 



i = l,..A. 



(29) 



The explicit expression for the denominator d and the coefficients hi are: 



d= {xl- xl) 
k 



2xiX2 ^ h x^ 7cosha;2t + ^ h xf^ X2 (2xi7coshxit + [xf + 7^) sinhxii) 

(30) 



Xl ( 1- ) (^2 + 7^) sinhx2t 



m 



X2 



X27 vo 



h X? — 2x9 ) + xo I h X? 17) coshx2t 

m / \m ' ' 



m 



+ I X2{vo — xo7)(7coshxit + xi sinhxit) 



+ ( vqxI ( h X2 

m 



k 



k 
m 



Xo ( — + x| ) X27 + Vo {-xl + X2) 7^ ) sinhx2t 



(31) 
(32) 
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&3,t = Xi[ h Xi 



— + Xo ) 7(-wo + xot) coshx2t 

TO / 

— — 2x? + Xo I + xo ( — + Xo ) 7 ) coshxi* 

TO / V TO / / 



— ^-t^oXi + xj^ + xqXi + xl^ 7 + fo (-Xi + xl) 7^^ sinhxii 



— xiX2 ( — + X? ) f 1- Xo 1 (fo — X07) sinhx2i 

TO / \ TO 



TO 



+ x? 



X1X2 ( -fo I — + X2 ) - Xo 



TO 



2xf + X2 ) 7 I coshxit 



+X2 (^—vo + x|^ 7 + Xo ^~7^ + ^1 (^1 - ^2 + '^^)^^ sinhxit 



+xi 1- X2 (wo — xo7)(x2 coshx2t + 7sinhx2t) 

. TO 



(33) 



(34) 
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